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ABSTRACT 

We analyze the dynamical expansion of the H II region, photodissociation 
region, and the swept-up shell, solving the UV- and FUV-radiative transfer, 
the thermal and chemical processes in the time-dependent hydrodynamics code. 
Following our previous paper, we investigate the time evolutions with various 
ambient number densities and central stars. Our calculations show that basic 
evolution is qualitatively similar among our models with different parameters. 
The molecular gas is finally accumulated in the shell, and the gravitational frag- 
mentation of the shell is generally expected. The quantitative differences among 
models are well understood with analytic scaling relations. The detailed physical 
and chemical structure of the shell is mainly determined by the incident FUV 
flux and the column density of the shell, which also follow the scaling relations. 
The time of shell-fragmentation, and the mass of the gathered molecular gas are 
sensitive to the ambient number density. In the case of a low density, the shell- 
fragmentation occurs over a longer timescale, and the accumulated molecular gas 
is more massive than in the case of a high density. The variations with different 
central stars are more moderate. The time of the shell-fragmentation differs by 
a factor of several with the various stars of M* = 12 — 101 M & . According to 
our numerical results, we conclude that the expanding H II region should be an 
efficient trigger for star formation in molecular clouds if the mass of the ambient 
molecular material is large enough. 
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1. Introduction 

Massive stars are ubiquitous in many star-forming regions, and their feedback processes 
have significant impacts on the surrounding interstellar medium (ISM). UV/FUV radiation 
photoionizes/photodissociates molecular gas, and the H II and photodissociation regions 
(PDR) expand around the massive star. The stellar winds and supernovae explosions give 
rise to expanding bubbles in ISM. The role of these feedback processes in relation to ISM 
has been studied from two competing perspectives on their effects on star-formation activity. 
The expanding H II region and bubble erode the parental molecular clouds and suppress star 
formation {negative feedback). On the other hand, the shock front (SF) emerges as the high- 
pressure region expands. The SF compresses the ISM and star formation can be triggered 
in very dense layers {positive feedback) . 

In this paper, we focus on the H II regions and their feedback processes. The role 
of the expanding H II regions is also discussed in terms of negative and positive feedback 
processes. Most of the massive stars are born in the giant molecular clouds (GMCs) and H II 
regions begin to expand in the dense cores and clumps in GMCs (ultra compact H II regions, 
e.g., Churchwell 2002). As the H II region expands, the molecular gas is gradually eroded. 
When the IF reaches the cloud surface, the ionized gas spouts from the cloud ("champagne 
flow", see Tenorio-Tagle 1979). The molecular material around massive stars is dispersed, 
so that star formation is quenched (e.g., Whitworth 1979; Franco et al. 1994). On the 
other hand, the SF emerges in front of the IF owing to the pressure excess of the H II 
region. The SF sweeps up the ISM to form a shell and star formation will be induced by 
the instability of the shell ("collect and collapse" scenario ; see, e.g., Elmegreen & Lada 
1977; Elmegreen 1989, and references therein). The SF also compresses pre-existing globules 
through their photoevaporation and star formation will be triggered (radiative implosion ; 
see, e.g., Bertoldi 1989; Bertoldi & McKee 1990). 

The hydrodynamical expansion of the H II regions has been studied numerically (e.g., 
Mathews 1965; Tenorio-Tagle 1979; Franco et al. 1990; Garcia-Segura & Franco 1998). These 
efforts have successfully shown the various dynamical aspects of the expansion in both ho- 
mogeneous and inhomogeneous ambient media. However, these do not include the outer 
PDR and thermal processes that should be dominant there. Hill & Hollenbach (1978) have 
analytically shown that the dissociation fronts (DFs) initially expand outside the SF, but 
finally slow down to be taken into the shell. Roger & Dewdney (1992) and Diaz-Miller et al. 
(1998) studied the expansion of the H II region and outer PDR around a single massive star, 
solving the radiative transfer of UV and FUV photons. Although the hydrodynamics has not 
been included, they have shown that a significant amount of the molecular gas is not ionized 
but photodissociated. Diaz-Miller et al. (1998) have pointed out that the photodissociation 
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of the molecular gas reduces the star-forming capacity of the cloud 

Deharveng et al. (2003) have provided a good simple example to investigate the role 
of the expanding H II region. They have observed the molecular fragments around the 
galactic H II region, Sharpless 104 (Shl04), and found a young cluster within the core of one 
fragment. They have argued that this is clear evidence of the "collect and collapse" model. In 
our previous paper (Hosokawa & Inutsuka 2005, hereafter, Paper I), we have modeled Shl04 
performing the hydrodynamical calculation including the PDR and the swept-up shell as well 
as the H II region. Our calculation has shown that the molecular shell forms around the H II 
region shielding the FUV radiation. The abundance of the accumulated molecules shows 
excellent agreement with the observation. We have also suggested a role for the expanding 
H II regions as a positive feedback process. The swept-up mass remains in the shell and the 
shell can be dominated by the cold molecular gas shielding UV and FUV photons. 

Recently, Deharveng et al. (2005) have reported dozens of galactic H II regions similar 
to Shl04, which show that such triggering is a common process in some situations. So, the 
next interesting questions to answer are the following : 

• How does the time evolution change with differing ambient number density and mass of 
the central star? Is the molecular gas always accumulated in the shell? If the evolution 
change with a different ambient density or central star, why so? 

To answer these questions, we investigate the time evolution of the H II region, PDR, and the 
swept-up shell with various ambient number densities and central stars. As in Paper I, the 
UV- and FUV-radiative transfer and the thermal and chemical processes are solved using a 
time-dependent hydrodynamics code. We show the numerical results of some representative 
models and explain what causes the differences among models. We can show that, in many 
cases, the FUV photons are shielded and the molecular gas is finally gathered in the shell. 

The organization of this paper is as follows: In §2, we explain the input physics of our 
code. The subsequent §3 is devoted to showing the procedures of the numerical calculation. 
§4 and 5 are the main part of this paper, where the numerical results are shown. The 
dependence on the ambient number density, and on the luminosity (mass) of the central star 
is studied in each section. §6 and 7 are assigned to discussions and conclusions. 
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2. Input Physics for Calculation 
2.1. Basic Equations 

£1.1. Hydrodynamic Equations 

In this paper, we consider the spherical expansion of the H II region. The hydrodynamic 
equations of the system are written as, 



dt \p J dm 

p&^ + ^W-a), (3) 



p = {0.5(1 + 3X H+ + X H ) + Vke} nfcT , (4) 

p = (1 + 4F Hc )^Hn = /im H n , (5) 

where m is the mass coordinate defined as dm = pr 2 dr, p is the mass density, n is the 
number density of the hydrogen nucleon, u is the velocity, p is the gas pressure, and E is the 
total energy, including the kinetic energy {u 2 /2) and the thermal energy (e =p/{^— l)/p) 
per unit mass. We adopt 7 = 5/3 for simplicity. In the above equations, Xh+ (Ah) is the 
number ratio defined as n H +/ n (%/n), and Y He is the abundance of He atoms. In the energy 
equation, T and A represent the heating and cooling rates per hydrogen nucleon, and q CO nd 
is the heat transfer by thermal conduction. We include 5 heating processes and 12 cooling 
processes to calculate T and A (see §2.2). The thermal conduction is needed to satisfy the 
"Field condition" (see Appendix B). We do not include the self-gravity of the gas to study 
the simple situation of the HII region expanding into a homogeneous medium. Even if we 
include the gas self-gravity, the H II region quickly expands at the velocity of ~ 10 km/s, and 
the smaller inward velocity caused by the self-gravity does not significantly affect our results. 
The gravity of the central star does not affect the evolution of H II regions, unless the H II 
region is very small. The radius where the velocity of escape from the star is comparable to 
the sound speed of the H II region is ~ 10~ 3 (M*/1OM ) pc, where M* is the stellar mass. 
This is much smaller than the initial Stromgren radius, 

with the number density considered in this paper, n < 10 4 cm -3 . 



The procedure to calculate a whole set of equations is explained in §3.1. 
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Table I. Properties of massive stars with solar metallicity (Diaz-Miller et al. 1998) 



T eff (10 4 K) a 


U. (M ) b 


log(S uv (s-l) c ) 


log(S FUV (s-l) d ) 


log(5 , Fuv/5'uv) 


5.0 


101.3 


49.89 


49.54 


-0.35 


4.0 


40.9 


48.78 


48.76 


-0.02 


3.3 


19.0 


47.75 


48.21 


0.54 


2.4 


11.7 


45.38 


47.39 


2.01 



a : effective temperature, b : stellar mass, c, d : emission rates of UV and FUV photons 



2.1.2. Radiative Transfer Equations 

We solve the radiative transfer of both UV and FUV photons. UV photons photoionize 
atoms in the H II region and heat up the gas by the photoionization. FUV photons pho- 
todissociate molecules in the PDR and heat up the gas mainly by the ionization of the dust 
grains (photoelectric heating). 

We solve the UV (hi/ >13.6 eV) and FUV (11.0 eV < hv < 13.6 eV) radiative transfer 
separately to avoid the time-consuming calculation. The UV radiative transfer equation is, 

1 <9(r 2 i^ uv ) 

r2 dr = ~ n ' ~ X h+)^,uv^,uv, (7) 

where F^uv du is the UV photon number flux at the frequency range, v — > v + du and a^uv 
is the cross section for the photoionization of hydrogen atoms. Equation (7) is solved as 

F v ,w = exp(-r I/)UV ), (8) 

where SV,uv dv is the photon number luminosity from the central star at a frequency range 
of v — > v + dv, and r^uv is the opacity written as, 

t„,uv = J n(l- X H +)cr l/)U v dr. (9) 

Diaz-Miller et al. (1998) provide the frequency-integrated photon number luminosities of 
massive stars in UV and FUV range. They have calculated with the line-bracketed LTE 
atmosphere models by Kurucz (1993). We use their results for the solar metallicity listed 
in Table 1. Note that some recent works present slightly larger values of Suv an d SVuv- 
For example, Suv given by Martins, Schaerer & Hiller (2005) is larger by a factor of ~ 
1.5 (3.0) for 20 M (40 M ) star than the adopted values. We have calculated the ratio, 
Suv/£fuv w ith the model spectrum distributed by Lanz & Hubeny (2003), and found that 
S F uy is also larger than the values in Table. 1 only by a factor of less than 2. However, 
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these differences by a factor of a few do not significantly affect our quantitative analysis. We 
approximate the spectrum of the Lyman-continuum, S u ^y, by the Planck function of the 
effective temperature, T eS and normalize S^uv so that the integrated UV luminosity agrees 
with the values listed in Table 1. We adopt the frequency-dependent cross section, 

<7„ )UV = 6.3 x 1(T 18 — (cm 2 ), (10) 

V^Lyl/ 

where z/L y i corresponds to the Lyman limit, hv\^\ = 13.6 eV. We do not take into account the 
absorption of ionizing photons by dust inside the HII region because of the uncertain size 
distribution of grains. We discuss the effect of the dust in the H II region in §6.1.1. The diffuse 
UV radiation produced by the recombination is treated by the on-the-spot approximation. 

The treatment of the FUV radiative transfer is a little more complicated. Hydrogen 
and carbon monoxide molecules are photodissociated by the line absorption of energetic 
photons (11 eV < hu < 13.6 eV). Lower energetic photons (hu < 11 eV) also contribute to 
some other processes. For example, photons of hu > 6 eV (typical work function of grains) 
contribute to the photoelectric heating and the dust recombination cooling (Bakes & Tielens 
1994). Below, we solve the transfer of three types of FUV flux separately, which are for the 
photodissociation of H 2 molecules, CO molecules, and for other microprocesses. 

We use the FUV flux, F v ^ 2 to calculate the photodissociation rate of hydrogen molecules 
with the self-shielding effect. The radiative transfer equation for -F„h 2 is, 

^ (r2 ^' H2) = -n(0.5X U2 a U:U2 + a d )F v ^ 2) (11) 

where ov,h 2 and era are the cross section of the photodissociation of hydrogen molecules and 
the dust absorption, Xh 2 is the molecular ratio defined as, Xh 2 = 2nn 2 /n = 1 — Xh+ — X-r. 
The dust attenuation law for the H 2 photodissociating photons is exp(— 2.5 Ay) (Tielens & 
Hollenbach 1985), and we adopt the corresponding cross section, a& = 1.2 x 10~ 21 cm 2 . 
Equation (11) is solved as, 

Fu,h 2 = ^^exp(-r i , i H 2 ), (12) 

where S u ^ 2 du is the FUV photon number luminosity at u — > u + du and r H2iJ , is the optical 
depth given by, 

t v ,h 2 = J n (O.5Xh 2 0V )H2 + <r d ) dr. (13) 

The photodissociation of H 2 molecules occurs via a so-called "Solomon process" (Stecher & 
Williams 1967). We assume that all H 2 molecules are at the rotational- vibrational levels 
of v — 0, J = (ortho) or v — 0, J = 1 (para) of the ground electronic level. The 
ortho-to-para ratio is assumed to be a typical value, 3:1, which is the limit for the high 



-7- 



temperature (Hollenbach & Salpeter 1970). We consider only transitions via Lyman-bands, 
which are the main channels for the "Solomon process" . In the Lyman-band transitions, the 
rotational quantum number changes as AJ = —1 (P transition) or AJ = +1 (R transition). 
Therefore, only 3 types of transitions are possible for each excited vibrational level in 5 X E+, 
v' ; (v , J) = (0,0) — > {y' , 1), (0,1) — > (f',0), and (0,1) — > (v',2). The vibrational quantum 
number, v and v', can change freely. For the line profile for each transition, we adopt a Voigt 
profile, 

7' f + °° ex P(-y 2 ) 

* M = FrfWiA^ * (14) 

where 7j is the natural line width, Az/d is the Doppler width, and v = [y — z/ )/Az/d, where 
u Q is the frequency of the line center. In our calculation, we use the constant temperature, 
T = 100 K to calculate the Voigt profile. We have confirmed that our numerical results are 
not sensitive to T or Az/d. As noted by Hollenbach et al. (1971), the column density of H2 
molecules easily exceeds 10 19 (cm -2 ), then the Lorentz wings dominate the absorption of 
FUV photons rather than Doppler cores (see Fig. 1 below). The frequency-dependent cross 
section of the line is, 

" M = jzifrM- (15) 

The total cross section, aj is given by <jj = ne 2 /m e c ■ fj, where fj is the oscillator strength. 
Numerical values of fj, 7j, and z/ are taken from Allison & Dalgarno (1970), Abgrall et 
al. (1992), and Abgrall et al. (1993). Although there are about 20 sets of three lines of 
Lyman-bands in 11.0 eV < hv < 13.6 eV, we approximate these many lines by adopting a 
representative set of lines. We pick up three v' = 12 lines. The logarithmic grids are set 
around these lines (12.82 eV < hv < 12.93 eV) so that the resolution becomes higher around 
each line center. We consider the constant radiation field, S v ^ 2 m this frequency range and 
normalize this so that the integrated luminosity is equal to the value listed in Table 1. The 
absorption cross section for Lyman-bands is calculated taking into account the line overlap, 



= (l-Zrem^ViM' ( 16 ) 



(see Fig. 1) where hj is the ratio of ortho or para H2 molecules and f rem is the probability 
of the reemission. We adopt / rem = 0.14 following Shull (1978). About another 15% of 
the pumped H2 molecules decay to be dissociated (see §2.1.3). We assume that the other 
pumped molecules decay to the excited (v, J) levels in the ground electronic state, followed 
by an infrared cascade. Fig. 1 shows the absorption cross section of the hydrogen molecules, 
ct Vj h 2 and the dust grains for which we adopted a&. As this figure shows, the effect of line 
overlap of Lyman-bands works at the column density, N H2 > 10 20 cm -2 . About half of FUV 
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photons are absorbed by dust grains rather than H 2 molecules, because o"d is larger than 
ov h 2 i n the feet of the Lorentz wings. 

The photodissociation of CO molecules is also caused by the line absorption of FUV 
photons (van Dishoeck & Black 1988). Since the distribution of such lines and the line 
overlapping with Lyman- Werner bands are complicated, we use other FUV flux, F^co to 
calculate the CO photodissociation rate with some shielding effects in the minimal frequency 
range. The radiative transfer equation for -F„ iCO is, 

~2 ~ ^ = ^ n ( X coZ c cru,co + 0.5Xh 2 (7h 2 ,co + cr dtC o)F UtC o, (17) 

where X co is the molecular fraction defined as, X co = rico/(Z c n), Z c is the abundance of 
C atoms, a^co, °"h 2 ,co, and o"d,co are cross sections for the CO photodissociation, shielding 
by H 2 molecules, and the dust absorption. Equation (17) is integrated as, 

K,co = exp(-7v,co), (18) 

where r^co is, 

r u ,co = J n {X co Z c a^co + 0.5X H2 o- H2 ,co + 0d,co) dr. (19) 

van Dishoeck & Black (1988) have listed about 30 lines to photodissociate CO molecules in 
the energy range of 11.0 eV < hv < 13.6 eV. As for the H 2 molecule Lyman bands, we use the 
typical one line, cr^co to calculate the photodissociation rate of CO molecules. We adopt the 
Voigt profile given by (14) and (15) again, and make up the representative line with typical 
parameters ; hv =12.925 eV (A = 960A), 7 = 1.0 x 10~ n s _1 and / = 0.01 respectively. 
We put the logarithmic grids at 12.88 eV < hv < 12.96 eV to get the resolution higher 
around the line center. The radiation field, S^co in equation (18) is determined in the same 
manner as S U: u 2 - 111 the PDR, CO molecules are shielded by the dust and H 2 molecules. Lee 
et al. (1996) have calculated the shielding function by the dust and H 2 molecules for the CO 
photodissociation rate. We have found that their results (Table 11) can be approximated 
with some analytic functions. Their dust attenuation law can be fitted with exp(— 3.5Ay) 
within a factor of 2 at Ay < 7. Similarly, their H 2 molecule shielding function can be fitted 
with exp(— 7Vh 2 /1-6 x 10 21 cm~ 2 ) within 20% error at A"h 2 < 8 x 10 21 cm" 2 . Therefore, 
we adopt the corresponding cross sections in our calculations, o" d CO = 1.7 x 10~ 21 cm 2 
and o"h 2 ,co = 6.25 x 10" 22 cm 2 respectively. Fig. 2 shows the cross sections for the CO 
photodissociation rate. As this figure shows, the self-shielding effect is not so significant 
owing to the small abundance of C atoms. The dust cross section is largest in most of the 
energy range, which means that CO molecules are mainly shielded by the dust absorption. 
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Finally, we define the frequency-integrated FUV photon number flux normalized with 
the Habing's unit, Gfuv (Habing 1968). We assume that Gfuv is attenuated only by the 
dust absorption, 

1 <9(r 2 G FUV ) 

r2 Or = _ncTdGFUV ' ( 20 ) 

which is solved as, 

GFUV = ^^^ eXp( " Td) ' (21) 
where F H is the Habing field, F H = 1.21 x 10 7 cm _2 s _1 (Habing 1968; Bertoldi & Draine 
1996), and r d is the dust optical depth defined as, 



= J na d dr. (22) 



We adopt cr d = 1.2 x 1CT 21 cm -2 , which is the average value between 1000 and 2000 A and 
corresponds to the attenuation law of exp(— 2.5Ay) (Roberge et al. 1991). As mentioned 
above, we assume cr d = in the H II region for simplicity. We use this flux, Gfuv to 
calculate the photoelectric heating and dust recombination cooling rate (Bakes & Tielens 
1994), reformation rate of CO molecules (Nelson & Langer 1997), and dust temperature 
(Hollenbach, Takahashi & Tielens 1991). 



2.1.3. Chemical Reaction Equations 

Non-equilibrium reaction equations are solved for the species of e, H, H + , H2, C + , and 
CO. The ionization rate of O is assumed to be the same as that of H. We adopt the solar 
elemental abundances with some depletions, F Hc = 0.1, Z Q — 5.4 x 10~ 4 , and Z c = 2.3 x 10~ 4 
(Hollenbach & McKee 1989). 

We solve the minimal set of chemical reactions listed in Table 2. The main reactions 
in the H II region are the photoionization by the UV photons from the central star and the 
recombination. In the PDR, the FUV photons photodissociate H 2 and CO molecules. The 
reformation of H 2 molecules mainly occurs on the surface of dust grains. The dissociation 
and reformation of CO molecules is approximated with the direct process between CO and 
C+ (Nelson & Langer 1997). 

Below, we show the photoionization and photodissociation rates, which relate to the 
transfer equations explained in §2.1.2. The photoionization rates of H and H 2 , photodisso- 
ciation rates of H 2 and CO molecules are, 

dt ~ df 



POD 

Xh / 0v,uv i^uv dv' , (23) 
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Table 2. Included Chemical Processes 



species 


reactions 


reference 


note 


H -> H+ 


photoionization 
collisional ionization 
cosmic-ray 


see §2.1.3 

1 

2 




H+ -> H 


recombination 


3 




H 2 -> H+ 


photoionization 
cosmic-ray 


see §2.1.3 
2 




H 2 -> H 


photodissociation 
cosmic-ray 


see §2.1.3 
2 




H -> H 2 


reformation 
gas-phase reaction 


4 
4 


n 


CO -> C+ 


photodissociation 


see §2.1.3 




C+ -> CO 


reformation 


5 





REFERENCES— (1) Tenorio-Tagle 1986, (2) Wolfire et al. 1995, (3) , (4) Hollenbach & 
McKee 1979 (HM79), (5) Nelson & Langer 1997 NOTES — fl : we use the rate by HM79 
for reactions; H + e — > H~ + 7 and H~ + H — > H 2 + e 



d ir = ~ 2 ^r = 2X112 [ ^' uv ^ uv ^' ' (24) 



rfXH -2^ = 2X Hi r + 7^-^^A/ ) (25) 



dt dt J v 1 — /] 



rft dt 



X C o I (?v,co F UtC o dv' , (26) 

J V- 



where pdis in equation (25) is the probability that the pumped H 2 molecules are photodisso- 
ciated, decaying to X 1 ^^. We adopt the average value, p^s — 0.15 (Draine & Bertoldi 1996). 
The integration in equations (25) and (26) are done in the considered frequency range around 
our representative lines. As equation (24) shows, we assume that if UV photons (> 13.6 eV) 
reach the molecular region, one H 2 molecule is directly ionized to two H + ions as two H 
atoms are photoionized. 



2.2. Thermal Processes 



We include 5 heating processes and 12 cooling processes in the energy equation. These 
processes are listed in Table 3. In Table 3, heating/cooling processes are divided into pro- 
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Table 3. Included Thermal Processes 





region 


process 


reference 


note 




HTT 


r4 nhnf ninTii7^finn 

LL IJLL(J HJUJlllZjdulwll 


1 pfr 








JT llULUClcCLlUll llUlil LILloL 


9 




Heating 


PDR 


H2 photodissociation 


3 








H 2 reformation 


3 








Cosmic-ray 


4 








H recombination 


1 etc. 








Lyman-a 


1 






HII 


01 (63.0/im) 


5 








Oil (37.29/im) 


5 








CII (23.26/im) 


5 




Cooling 




Collisional ionization 


6 








01 (63.1 jum) 


5, 7 


n 






CII (157.7/im) 


5,7 


n 




PDR 


H 2 rot/vib excitation 


3, 8 








CO rot/vib excitation 


3 








Dust recombination 


2 








Collisional dust-gas heat transfer 


5 


t2 



REFERENCES— (1) Spitzer 1978, (2) Bakes & Tielens 1994, (3) Hollenbach & McKee 1977, 
(4) Shull & Van Steenberg 1985, (5) Hollenbach & McKee 1989, (6) Tenorio-Tagle 1986, (7) 
Tielens & Hollenbach 1985, (8) Galli & Palla 1998. NOTES — fl : these can be a heating 
process when the background radiation is greater than the source function, f 2 : this can be 
a heating process when the dust temperature is higher than the gas temperature. 
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cesses dominant in the H II region or PDR. In the H II region, the included heating process 
is photoionization heating of hydrogen atoms. The main cooling processes in the H II region 
are recombination cooling and some meta-stable atomic lines (e.g., [O II] 3729 A, [C II] 2326 
A). We do not include the photoelectric heating in the H II region, because this is sensitive 
to the uncertain abundance of small grains. We do not include the line cooling with highly 
ionized ions (e.g., [O III], [C III] etc.) either. We are interested in the boundary region 
between the H II region and the molecular cloud, and thus we do not solve the detailed 
ionization structure in the H II region. 

In the PDR, the main heating process is photoelectric heating. We separately list the 
dust recombination cooling, which exceeds the photoelectric heating at high temperatures 
(T > 10 3 K, Bakes & Tielens 1994). In our calculation, cosmic-ray heating is important 
only in the outer molecular region, where the FUV flux from the central star is significantly 
reduced by geometrical dilution and dust extinction. The main cooling processes in the 
PDR are fine-structure atomic lines of [O I] 63.1/xm and [C II] 157.7/im, which are the 
characteristic lines of the PDR (e.g., Hollenbach & Tielens 1999). The [O I] 63.1/im is more 
efficient than [C II] 157.7/im at high density (n > 10 4 cm" 3 ), and it is the reverse at the 
low density (n < 10 2 cm" 3 ). We do not include [CI] 370/im and 609/zm transitions because 
we adopt the approximate treatment that CO molecules are directly photodissociated to C + 
(Nelson & Langer 1997). The rotational transitions of CO molecules are the primary cooling 
process in the molecular region. Collisional heat transfer between gas and dust is important 
when the density is very high (> 10 5 cm -3 ). If dust temperature is higher than the gas 
temperature, this process heats up the gas. 



2.3. Dust Temperature 

We calculate the dust temperature following the approximation method by Hollenbach, 
Takahashi & Tielens (1991). In our calculation, several thermal processes depend on the 
dust temperature. Dust IR emission affects the level population of atoms and cooling rates 
with line transitions (e.g., [01] 63.1/im and [C II] 157.7/mi). If the intensity of the IR dust 
emission at the line frequency exceeds the source function of the line, this transition works 
as a heating process rather than a cooling process. The collisional dust-gas heat transfer 
rate is proportional to the difference between the gas and dust temperature. 
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3. Calculation Procedures 
3.1. Outline 

Our numerical scheme for the hydrodynamics is based on the 2nd-order Lagrangian 
Godunov method (see, e.g., van Leer 1979). First, we explicitly integrate the continuity 
equation and the momentum equation to obtain the new position of the cell, r - n+1 , gas 
density, p™ +1 , and velocity, Uj n+1 , where the indices j and n refer to different steps in space 
and time respectively. Next, we iteratively solve the chemical rate equations and the energy 
equation. The energy equation (3) is transformed with equations (2), (4) and (5) as, 



dT T /dX u 3 dX n+ \ 

dt 2w V dt dt J 



/im H (7 ~ 1) f 2 d P d(r 2 up) d(r 2 q cond ) \ 
kw \ dm dm dm J 

= ^<r-A). (27) 



where w is defined as, 



w = 0.5(1 + X H + 3X K+ ) + Y He . (28) 



The terms with m-derivative are explicitly calculated before the iteration. The different rate 
equations and energy equation include the new gas temperature, Tj n+1 , dust temperature, 
Tdj ra+1 , and chemical compositions, X s j n+1 (s : species). We determine these quantities 
by means of the Newton-Raphson method. After the iteration converges, the gas pressure 
and UV/FUV flux are calculated using the newly obtained gas density, temperature, and 
compositions. The integration proceeds from the center to the outer cells. 



3.2. Initial Condition 

The hydrodynamical expansion of the H II region occurs after the IF passes the initial 
Stromgren radius (expansion phase). We start our calculation when the H II region is still 
in the previous formation phase. In the formation phase, the timescale of the system is very 
short and the expanding H II region hardly affects the gas dynamics. At t = 0, we set the 
constant mass density and no initial velocity field. We set the initial H II region (A H + = 
X c+ = 1, T = 10 4 K) at r < R st /5. The outer initial PDR (X H = X c + = 1, T = 100 K) 
and the molecular region (X H2 = X C o = 1, T = 10 K) is set at R st /5 < r < 2R st /5 and 
r > 2R st /5 respectively. 
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Table 4. Model Parameters 



Model 


M* (M ) a 


n H ,o(cm" 3 ) fe 


Rst (pc) c 


t dyn (Myi) d 


S101 


101.3 


10 3 


1.27 


0.12 


S41 e 


40.9 


10 3 


0.56 


0.05 


LD-S41 


40.9 


10 2 


2.53 


0.23 


HD-S41 


40.9 


10 4 


0.12 


0.011 


S19 


19.0 


10 3 


0.25 


0.023 


S12 


11.7 


10 3 


0.04 


0.004 



a : mass of the central star, b : number density of the ambient gas, c : initial Stromgren 
radius, d : dynamical time, R st /Cmi, where Ccn is the sound speed at T = 10 4 K. e : Our 
model for Shl04 (Paper I). 



3.3. Time Step 

The time step of the calculation, At, is constrained by some conditions. The first one 
is the normal Courant condition, 

My = p- , (29) 

where C s j is the sound speed at the j-th cell. We include the thermal conduction explicitly 
in the energy equation, and we have another Courant condition, 



^ = kuMTT)^ (30) 



where Kj is the conductivity calculated with equation (B2) to satisfy the Field condition. 
Besides these Courant conditions, At is constrained by the velocity of the IF. In the early 
formation phase of the H II region, the IF expands supersonically. Since the ionization rate 
and the temperature change sharply at the IF, At must be short enough that the IF advances 
less than the width of one cell, Arj in one time step, At. 

We choose the minimum time step from among the above time steps to determine At. 
In the early formation phase, At is mainly constrained by the expanding velocity of the IF. 
In the next expansion phase, when the SF emerges, the Courant conditions limit the time 
step. 
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4. Dependence on Ambient Number Density 
4.1. Models 

In this section, we show how the time evolution of the H II region, PDR, and the swept- 
up shell changes with different ambient number densities. For comparison, we calculate 
the expansion of the H II region around the 41M Q star, which is the central star of Shl04 
modeled in Paper I, but with different ambient number densities. As we have adopted the 
ambient number density of n^o — 10 3 cm -3 in Paper I, we show the results of the higher 
density case, nn,o — 10 4 cm -3 and the lower density case, nn,o — 10 2 cm -3 below. The initial 
Stromgren radius and the dynamical time for the calculated models are listed in Table 4. 
All the other conditions are the same as those in our model of Shl04 (S41 in Table 4) except 
for the number density. 



4.2. Scaling Relations 

First, we show some scaling relations that are very useful to interpret the numerical 
results. The initial Stromgren radius depends on the ambient number density, and other 
relevant quantities are also scaled as, 

Rst oc n- 2 / 3 , (31) 

*dyn = Rst/C mi oc n- 2 / 3 , (32) 

M sh oc nR 3 st oc n" 1 , (33) 

N sh oc nR st oc n 1/3 , (34) 

F FUV ,i oc S FVV /Rl oc n 4 / 3 , (35) 

where M s h and is the mass and the column density of the shell, -Fpuv.i is the incident 
FUV flux at the IF at a given t/t dyD . To derive the above relations, we assume that most 
of the swept-up gas remains in the shell, which is a good approximation for expansion in 
a homogeneous ambient medium. As (31) and (32) show, the larger H II region expands 
in the lower density ambient medium over a longer timescale. Equation (33) shows that 
the swept-up mass is larger with the lower density medium. The shell structure is mainly 
determined by the column density of the shell, N s ^ and the FUV flux at the IF, -Fpuv.i- The 
detailed properties of the shell are explained with equation (34) and (35) in §4.3. 
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4.3. Results of Numerical Calculations 

Figures 3 and 4 show the hydrodynamical evolution of models HD-S41 and LD-S41 
respectively. Since the Stromgren radius is smaller for the higher ambient density, the 
H II region expands to the region of ~ 1 pc in HD-S41, and ~ 20 pc in model LD-S41 in 
t ~ 20 t dyn . As expected with equation (32), the timescale of the evolution is longer for 
the lower density model. The time of 20 t dyn corresponds to ~ 0.25 Myr for HD-S41, and 
~ 5 Myr for model LD-S41. 

The basic time evolution is the same as that of model S41 (Paper I). The SF emerges 
when the IF reaches the initial Stromgren radius. The gas density within the shell is about 
10-100 times as high as that of the ambient medium, though the detailed structure of the 
shell differs between HD-S41 and LD-S41 (explained below). The time evolution after the 
SF emerges is also similar to that of model S41. The time evolution of the radial expansion 
is well approximated with, 

'W'H^ + lVf^)'. (36) 
which is derived from the equation of motion of the shell, 

j t (f J&Poite) = ^Rlmi = ^RhPoC^ , (37) 

where po is the mass density of the ambient medium, and pn is the gas pressure in the 
H II region. Equation (36) is different from the well-known expansion law (e.g., Spitzer 
1978) by a factor of a/4/3, but shows better agreement with the numerical results. As 
equation (37) shows, equation (36) does not involve the detailed structure and the chemical 
composition of the shell, which we analyze with our numerical calculations. The swept-up 
mass in the shell is larger for the lower ambient number density, as equation (33) shows. For 
example, M sh ~ 1O 5 M in model LD-S41 and M sh ~ 1O 3 M in model HD-S41 at t ~ 20 t dyn 
respectively. The upper panels of Fig. 5 and 6 show the time evolution of various front 
positions. Initially, the DFs expand more rapidly than the SF and the IF, then the PDR 
appears in front of the SF. As the H II region expands, the IF and the SF gradually overtake 
the preceding DFs of H 2 and CO in turn. Finally the DFs are taken up into the shell. The 
SF enters the molecular region, and the molecular gas is accumulated within the shell. The 
rapid reformation within the shell also helps to accelerate the accumulation of the molecular 
gas. 

The main difference among these models is the time (or t/t dyn ) when the DFs are taken 
into the shell. In the higher density model, HD-S41, the IF and the SF catch up with both 
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DFs of H 2 and CO by the time of 5 td yn - In the lower density model, LD-S41, the DF of 
H 2 is quickly taken into the shell, but the accumulation of CO molecules is delayed. The 
DF of CO is taken into the shell at t ~ 20 - 30 t dyn in model LD-S41. About 70% (90%) of 
carbon (hydrogen) atoms within the shell exist as CO (H 2 ) molecules in model HD-S41 at 
t ~ 20 tdyru but 20% (90%) in model LD-S41 at the same time step. We explain the reasons 
for this difference below. With the higher ambient number density, the incident FUV flux 
at the IF is strong, as equation (35) shows. Although the FUV photons are available to 
photodissociate molecules, our calculation shows that the FUV photons are quickly shielded 
in model HD-S41. This is because the column density of the shell, N sh is larger for the higher 
density model, as equation (34) shows. Fig. 7 is a schematic picture showing the density- 
dependence of Agh and -F F uv,i- The incident FUV flux is attenuated enough by the dust 
extinction through the shell, and the molecular gas is reliably protected in the outer part of 
the shell. Though the density-dependence of A s h, A s h oc n 1 ^ 3 , is weaker than that of -F F uv,i> 
^Fuv,i ^ n 4 / 3 , the dust absorption becomes significant with a high number density owing 
to the attenuation law of exp(— <7dA sn ). The lower panels of Fig.5 shows the time evolution 
of the column density of each region in model HD-S41. The column density of the shell at 
t ~ 5 t dy n is already A sh ~4x 10 21 cm~ 2 (A v , s h ~ 2) in this model. In model LD-S41, on the 
other hand, the incident FUV flux is weak, but the column density of the shell is low. Since 
the CO molecules are mainly shielded by the dust absorption, it is not until A sn becomes 
somewhat large (Av, s h ~ 1) that CO molecules are accumulated in the shell. However, even 
in the early phase when Av jB h < 1, H 2 molecules can be quickly accumulated in the shell, as 
Fig. 6 shows. This is due to the efficient self-shielding of H 2 molecules. Hollenbach & Tielens 
(1999) have briefly shown that H 2 molecules can be self-shielded, if G-p\jv/ n < 4 x 10~ 2 
at the IF. In our calculations, Gpuv/n ~ 1-6 x 10~ 2 in model LD-S41 and Gfuv/™ ~ 0.5 
in model HD-S41 at t ~ 5 tdyn- Therefore, H 2 molecules are mainly self-shielded in model 
LD-S41, and shielded by the dust in model HD-S41 respectively. The self-shielding effect for 
CO molecules is not so efficient, owing to the small abundance of carbon atoms. 

The detailed structure of the shell differs significantly among these models. Fig. 8 (Fig.9) 
shows the detailed structure of the density and the temperature within the shell at t ~ 15 td yn 
in model HD-S41 (LD-S41). These figures show that the gas temperature in the inner region 
of the shell is higher in the model with the higher ambient density; Tfii ~ 500 K in model 
HD-S41 and T ffl ~ 50 K in model LD-S41. This is because the FUV flux at the IF, F FU v,i 
is larger for the higher ambient density owing to the smaller Stromgren radius (see eq.(35)). 
The FUV radiation field, Gfuv defined by eq.(21) at the IF is Gfuv ~ 10 in LD-S41 and 
~ 5000 in HD-S41 at the time step of both figures, t ~ 15 td yn - Within the shell, the 
temperature has a negative gradient and the density has a positive gradient, which is similar 
to model S41 (Paper I). Figs. 8 and 9 show that the changes in gas temperature and density 
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are milder in model LD-S41. The gas temperature decreases from 500 K to 25 K through 
the shell in model HD-S41 and from 50 K to 25 K in model LD-S41. These differences are 
also explained in equation (34) and (35). In model HD-S41, the FUV flux is strong and the 
column density is high. The intense FUV radiation is shielded by the shell with the large 
column density. In model LD-S41, the incident FUV flux at the IF is relatively weak and 
the column density of the shell is low. The gas temperature within the shell is lower and 
nearly constant. The FUV flux is attenuated from Gfuv ~ 5000 to 0.1 through the shell in 
HD-S41, and from G FUV ~ 10 to 1.5 in LD-S41 at the time of Fig.8 and 9. In model HD-S41, 
the gas temperature significantly decreases toward the outer part of the shell, because the 
photoelectric heating becomes inefficient. Conversely, the gas density of the outer region of 
the shell rises to n ~ 10 6 cm" 3 . The collisional heat transfer between the gas and dust is 
efficient there, and the gas temperature is almost equal to the dust temperature. The the 
swept-up shell is geometrically thicker in model HD-S41. Figs. 8 and 9 also show that the 
density jump at the IF is larger with the lower ambient density. This can be explained by 
the analytic jump condition at the weak D-type IF, 



Phii C 1 Ti 



1 + 5) (38) 



Phi 2C n T] 



HII 



(Mihalas & Mihalas 1984), where C\ and Cn are the sound speed ahead and behind the 
IF, and 5 is the quantity less than unity. Since the temperature T H i is lower for the lower 
ambient density, the density at the inner edge of the shell is higher in model LD-S41. 



4.4. Fragmentation of the Shell 

It is important to know whether or not the swept-up molecular layer fragments, as the 
triggered star formation scenario predicts. We do not consider the dynamical instability, but 
only the gravitational instability for the fragmentation of the shell, as discussed in Paper I. 
In the lower panels of Figs. 5 and 6, the shaded region is the expected unstable region, where 
t > (Gp)" 1 / 2 . The swept-up shell in each model is expected to suffer from gravitational 
instability. The unstable region spreads from the outer edge of the shell, which is the same 
as in model S41 (Paper I). In model HD-S41, the density is relatively low in the inner part 
of the shell. The unstable region appears at t ~ 0.08 Myr (7.3 tdyn) and gradually extends 
from the cold molecular layer, as Fig. 5 shows. As explained in the previous subsection, both 
H 2 and CO molecules are shielded by the dust, and the shell is dominated by these molecules 
when the shell becomes unstable. In model LD-S41, the swept-up shell is relatively cold and 
dense (see Fig. 9). The unstable region rapidly spreads in the whole shell at t ~ 1.5 Myr 
(6.5 tdyn)- In this model, the shell is dominated by H 2 molecules due to the efficient self- 
shielding effect by the time of the shell-fragmentation. However, the column density of the 
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shell is still not high enough to protect CO molecules by the dust extinction in the early 
expansion phase. The accumulation of CO molecules begins after the shell becomes unstable, 
which is different from model HD-S41. Regardless of some differences, the expanding H II 
region can sweep up the molecular gas, and the fragmentation of the shell is finally expected 
with different number densities. The star formation will be triggered in the fragmented shell 
along the lines of a "collect and collapse" scenario. 

5. Dependence on the Central Star 
5.1. Models 

In this section, we investigate how the time evolution of the H II region, PDR, and 
the shell depends on the UV/FUV luminosity (mass) of the central star. For this purpose, 
we calculate the expansion around the central star of 101. 3M , 19.OM and 11.7M as 
representative cases. We use the same ambient number density, nn,o ~ 10 3 cm -3 , as in 
model S41 (Paper I). Model parameters of these calculated models are listed in Table 4. 



5.2. Scaling Relations 

The less massive central star has smaller UV and FUV luminosities, but a large lumi- 
nosity ratio, Sfuv/Suv (see Table 1). The Stromgren radius depends on the UV photon 
number luminosity, Sjjy, as well as the number density, so that other relevant quantities also 
depend on Suv- The scaling relations which are similar to equation (31) - (35) are derived 
as, 

Rst oc S 1 ^ , (39) 

1 /3 

tdyn = Rst/ Cmi OC Sjjy , (40) 

M sh oc nR 3 st oc S uv , (41) 

N sh oc nR st oc , (42) 

oc S Fvy /R 2 st oc ( W/Suv) 2/3 Sfuv, (43) 

where M sh , N sh , and -F FU v,i are the quantities at a given t/td yn - As equation (39) and (40) 
show, the H II region expands into the larger region on a longer timescale, as the mass of 
the central star increases. The FUV flux at the IF, Fpuv,i depends on both the luminosity 
ratio, Sfuv/*Suv an d the FUV luminosity, Spuv- Since the dependence on the luminosity 
ratio is stronger than that on the FUV luminosity, the FUV flux, -Ffuv.i becomes larger for 
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the lower mass central star. The calculated shell structure reflects the column density, N sh 
(eq.(42)) and the FUV flux at the IF, F FUV ,i, (eq.(43)) which is explained in §5.3. 

5.3. Results of the Numerical Calculation 

Figs. 10, 11, and 12 show the hydrodynamical evolution of models S101, S19, and S12 
respectively. The basic time evolution features are similar among these models. When the 
IF reaches the initial Stromgren radius, the SF emerges in front of the IF. In all these 
models, the H II region expands into the region of ~ 1 — 10 pc over a timescale of ~ 1 Myr. 
However, the time normalized with td yn differs considerably among these models, because 
tdyn decreases as the UV luminosity decreases (eq. (40)). For example, the time of 1 Myr 
corresponds to ~ 8 tdyn in model S101, ~ 40 tdyn in model S19, and ~ 300 tdyn in model S12. 
Figs. 10 - 12 all show that the PDR initially appears in front of the SF and is finally taken 
into the shell over a timescale of ~ 1 — 2 Myr; such evolution is largely different in terms of 
time normalized with t dy n- 

These differences are explained using equations (42) and (43). As equation (43) shows, 
the FUV flux at the IF at a given t/tdyn, ^Fuv,i is larger with the less massive star. Conversely, 
the column density of the shell, JV B h is smaller for the less massive star, shown in equation 
(42). At t ~ 10 t dyn , for example, G FUV ~ 100 at the IF and A" sh ~ 3 x 10 21 cm" 2 in 
model S19, and G FVY ~ 600 and N sh ~ 7 x 10 20 cm~ 2 in model S12. Therefore, the 
strong FUV radiation is not shielded, and the reformed molecules are easily photodissociated. 
Fig. 16 is the schematic picture showing the luminosity-dependence of N s ^ and Fpuv.i- The 
accumulation of molecular gas in the shell is significantly delayed on t/t dyn . The upper panels 
of Figs. 13 - 15 show the time evolution of the various front positions in each model. Both 
DFs of H-2 and CO molecules are taken into the shell by the time of 4 td yn in model S101, 
45 tdyn in model S19, and 550 td yn in model S12. The times of the front overtaking are totally 
different at t/t dyn , but these times differ only by a factor of a few at time t; t ~ 0.4 Myr 
in model S101, t ~ 1 Myr in model S19, and ~ 2 Myr in model S12. Eventually, sufficient 
molecular gas always gathers in the outer region of the shell in a few Myrs with the ambient 
number density of nn,o ~ 10 3 cm -3 . The total mass of the swept-up shell (between the IF 
and the SF) is 3.5 x 10 4 M , 9.8 x 10 3 M , and 3.8 x 10 3 M© in models S101, S19, and 
S12 respectively at the final time step of Figs. 13 - 15. At that time step, more than 85% of 
hydrogen atoms are included in H 2 molecules in all models. The fraction of CO molecules is 
55%, 40%, and 20% in each respective model. 

The upper panel of Fig. 17 shows the time evolution of the maximum density of the shell 
in models S101, S41, S19 and S12. As this figure shows, the swept-up shell is always denser 
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with the more massive central star. This is explained with the time evolution of the velocity 
of the shell. The expansion of the H II region approximately follows equation (36), and the 
expansion velocity decreases as, 



because the gas pressure of the H II region decreases. In model S101, the FUV radiation is 
easily shielded because of the high column density of the shell. The photoelectric heating 
outside the shell becomes inefficient, and the gas temperature in front of the SF quickly 
decreases. At that time, t/tdyn is still small, and the expansion velocity is still large according 
to equation (44). The Mach number of the SF is large, Ai ~ 10. In model S12, it takes 
hundreds of tdyn to shield the FUV radiation by the dust extinction through the shell. The 
expansion velocity significantly decreases by (44), and the Mach number is only about a few. 
This difference causes a different shell density in each model, because the density jump at the 
isothermal SF is proportional to M 2 . The lower panel of Fig. 17 presents the time evolution 
of the Mach number in each model. This figure clearly shows that the Mach number is larger 
in the model with the more massive central star. 

During the expansion, the pressure gradient in the PDR outside the shell hardly affects 
the hydrodynamics. The pressure gradient outside the shell gradually disappears as the FUV 
photons are shielded by the dust in the shell. The SF does not emerge in front of the DFs 
in our calculations. 



Even with different central stars, the molecular gas is finally accumulated in the swept-up 
shell. The lower panels of Figs. 13, 14, and 15 show the time evolution of the column density 
of each region in models S101, S19, and S12 respectively. In all these models, the number 
density of the molecular layer finally becomes over ten times as high as the ambient density, 
and gravitational fragmentation is expected there. In Figs. 13 - 15, the shaded region is the 
unstable region expected in each model, where t > (Gp)^ 1 ^ 2 . As shown in these figures, the 
unstable region spreads in the shell, as the DFs of CO are taken into the shell. The unstable 
region appears at t ~ 0.2 Myr in model S101, t ~ 0.5 Myr in model S19, and t ~ 1.1 Myr 
in model S12. As explained in the previous subsections, these times are largely different at 
t/t dyn ; ~ 2 in model S101, ~ 20 in model S19, and ~ 300 in model S12. The FUV flux at 
the IF is stronger with the less massive star, and it takes a long time in t^yn units to form 
the dense molecular layer shielding the FUV radiation field. 




(44) 



5.4. Fragmentation of the Shell 
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6. Discussions 
6.1. Dust Grains 

6.1.1. Dust Grains in H II Region 

In the above calculations, we have not included the dust grains in the H II regions 
to avoid uncertainties (e.g., size distribution), but some observations indicate their presence 
(e.g., see Spitzer 1978). The dust absorption of the UV and FUV radiation will decelerate the 
propagation of the IF and DFs and affect the structure of the swept-up shell. Weingartner 
& Draine (2001) have noted that the photoelectric heating can be as important as the 
photoionization heating, even in the H II region. The importance of the dust absorption in 
the H II region is often evaluated with the column density of the dust-free static H II region, 

N mi oc nR st oc n 1/3 4v (45) 

(Petrosian et al. 1972; Arthur et al. 2004), where R st is the dust-free Stromgren radius, n is 
the ambient number density, and Suv is the UV photon number luminosity of the central 
star. The radius of the corresponding dusty H II region is smaller for larger Arii- Arthur 
et al. (2004) has also shown that the dust absorption in the H II region gradually becomes 
inefficient, as the H II region dynamically expands. 

In order to evaluate the effect of the dust in the H II region, we have calculated some 
models including the dust in the H II region as well as the outer PDR. In models S41-HD 
and S41-LD, the radius of the H II region and PDR decreases only by 10% and 3% compared 
with the dust-free counterparts. In models S101 and S41, the reduction of each region is 20% 
and 5% of the dust-free radius. These results agree with equation (45), which shows that 
the dust absorption in the H II region becomes significant with the higher ambient number 
density, or the more massive central star. In all of the above calculations, however, the dust 
absorption in the H II region is not so significant. Furthermore, the dust in the H II region 
also absorbs FUV photons in addition to UV photons. As briefly mentioned in Paper I, this 
only slightly promotes the accumulation of molecules in the shell. Therefore, dust grains in 
the H II region hardly affect our conclusions. The dust absorption in the H II region is more 
important in the denser and more compact region (e.g., ultra compact H II region), where 
the detailed treatment of grains is required. 
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6.1.2. Reformation of Hydrogen Molecules on Grain Surface 

The reformation of H2 molecules on the grain surface is one of the key processes in our 
calculations. The FUV radiation is consumed in photodissociating the reformed molecules 
in the shell, which accelerates the accumulation of H 2 molecules. However, there still remain 
uncertainties in the reformation rate, which depends on the detailed behavior of H atoms 
and H 2 molecules on the grain surface. At high dust temperatures, for example, H atoms 
absorbed on the grain surface evaporate before the reformation of H 2 molecules. The critical 
temperature, T cr , is about 20 K on pure grains, but Hollenbach & Salpeter (1971) have shown 
that the irregular surface with impurity sites raises T cr to 25 K < T cr < 50 K. Recently, 
Cazaux & Tielens (2002, 2004) have refined the absorption process of H atoms on the grain 
surface, and obtained a much higher critical temperature of T cr ~ 100 K. In our calculations 
presented in §4 and 5, we have not included the high-temperature cut-off in the reformation 
rate, for simplicity. We analyze the validity of this treatment here. 

We calculate some models including the cut-off using the functional form given by 
Hollenbach & McKee (1979) (HM79). In each model, two cases of T cr = 65 K (HM79) and 
35 K are applied. With T cr = 65 K, however, the H 2 molecular abundance in the shell is 
hardly affected in all examined models. With T cr = 35 K, the accumulation of H 2 molecules 
in the shell is slightly delayed, but the shell is finally dominated by H 2 molecules. The H 2 
abundance in the shell is reduced only by 2% in model S41 (Paper I) and 5% in model S41- 
HD at t ~ 20 tdyn- In the region where H 2 molecules accumulate, the FUV radiation is well 
shielded and the dust temperature is low. Therefore, even if the reformation is inefficient 
above several ten K, our conclusions do not change. 

6.1.3. Small Grains 

Our calculations do not include detailed energy and ionization balance of the small grains 
such as polycyclic aromatic hydrocarbons (PAHs). The ionization rate and size distribution 
of small grains are important factors for the PDR structure. The photoelectric heating rate 
depends on the abundance and ionization rate of PAHs, and the dust attenuation law of 
FUV photons may be slightly affected by the change of the grain size distribution. We have 
calculated some models with different photoelectric heating rates given by Weingartner & 
Draine (2001) based on the different abundance of small grains. The temperature profile 
in the PDR would be modified by a factor of less than 2, but this does not affect our 
conclusions. Some recent observations detect the strong PAH emission around the H II 
region. This emission sometimes shows the clear ring-like structure, which should trace the 
shell (e.g., Zavagno et al. 2006). The expected PAH emission related to the time evolution 
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of the shell will be very useful for these observations, which should be included in future 
modeling. 

6.2. Realistic Structure of Molecular Clouds 

The actual molecular clouds have complex structures rather than the homogeneous 
density distribution assumed in this paper. Below, we briefly discuss how realistic density 
structures affect the time evolution of the swept-up shell, PDR, as well as the H II region. 

6.2.1. Density Gradient 

When the H II region forms near the edge of a molecular cloud, or in a cloud with 
the density gradient (e.g., n oc r~ w ), its time evolution is different from that in a homoge- 
neous medium. The high-pressure ionized gas is hardly trapped, and flows out of the cloud 
(Tenorio-Tagle 1979; Franco et al. 1990). Such "blister-like" or "champagne flow" features 
have been observed in some H II regions. In these regions, only a small fraction of the swept- 
up gas remains in the shell, and the expanding IF finally destroys the parental molecular 
cloud. As argued by Whitworth (1979) and Franco et al. (1994), the rapid ionization of 
the molecular gas will quench the star formation activity in the clouds. Furthermore, we 
note that the star-forming capacity of the clouds will be limited by FUV photons as well as 
UV photons. Unless the swept-up gas remains in the shell, stellar FUV radiation is hardly 
attenuated through the shell. Such unshielded FUV photons dissociate molecules efficiently 
over a much wider region than the H II region, which reduces the star formation efficiency of 
the clouds. The dynamical expansion of the H II region, PDR, and the shell in a cloud with 
a density gradient has not been studied well, and this will be explored in our forthcoming 
papers. 

6.2.2. dumpiness 

Real molecular clouds generally have clumpy or hierarchical structures. The pre-existing 
density inhomogeneity in molecular clouds can modify the structure and the evolution of the 
H II region, PDR, and the swept-up shell. 

In the H II region, the pre-existing clumps are gradually photoevaporated by stellar 
UV radiation. These clumps are basically smoothed out in the sound crossing time of each 
clump, which is much shorter than the dynamical time of the H II region, td yn = -Rst/Cmi- 
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Outside the H II region, the PDR structure depends on how deep the FUV radiation 
penetrates into the molecular region. Some previous work (e.g., Boisse 1990) has shown 
that the radiative transfer in a clumpy medium can be reduced to that of a homogeneous 
medium with the same average density and lower effective extinction. In the PDR around 
the H II region, the FUV photons will easily seep into the molecular region, threading the 
clumpy structure. This will delay the SF and IF overtaking the DFs, and the accumulation 
of molecules in the shell. In clumpy molecular clouds, on the other hand, most of the gas is 
thought to be contained in optically thick clumps. Molecules in these opaque clumps will be 
protected against the FUV radiation by the self-shielding of these clumps. In order to study 
these two competing effects, mult i- dimensional and time-dependent calculations are needed. 
Note, however, that our 1-D calculation has successfully explained the observed molecular 
abundance around Shl04 (Paper I). Thus we expect that our 1-D calculations should be a 
good approximation for the averaged physical quantities. 



6.2.3. Turbulent Velocity Field 

Observation shows that there are the supersonic turbulent motions in molecular clouds 
(e.g., Elmegreen & Scalo 2004). Since the velocity dispersion increases with the lengthscale 
(~ a few km/s at about 10 pc), the effect becomes important with the lower ambient density. 
The expanding shell actually dissolves when the shell velocity becomes less than the velocity 
dispersion of the turbulent motion. In §5, we have shown that the shell velocity becomes 
less than a few km/s when the molecular gas accumulates in the shell with the lower-mass 
star. In this case, the shell will dissolve without gathering the molecular material, and the 
triggered star formation will not occur. We study the net feedback effect in the turbulent 
molecular cloud in a subsequent paper. 



6.3. After the Cloud Destruction 

In this paper, we have considered only the expanding H II region in the molecular 
cloud, which can trigger the propagation of star formation activity within the cloud. After 
the molecular cloud disperses due to the negative feedback from massive stars and/or the 
external turbulent motion, stellar UV photons spread into the diffuse neutral medium. The 
H II regions expand in the ambient neutral medium, but the SF preceding the IF should 
also compress the neutral gas to enable the rapid reformation of molecules (e.g., Koyama 
& Inutsuka 2000; Bergin 2004, and references therein). The role of the giant H II regions 
expanding in the diffuse neutral medium will be explored in our next paper, Hosokawa & 
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Inutsuka (2006, in preparation). 

7. Conclusions 

In this paper, we have studied the time evolution of the H II region, PDR, and the 
swept-up shell around a massive star, performing the numerical calculation of radiation- 
hydrodynamics. Following Paper I, we have investigated how the time evolution changes 
with the different ambient number density and the luminosity (mass) of the central star. The 
basic evolution with different parameters is qualitatively similar to the fiducial one explained 
in Paper I, but quantitatively different. We have derived some simple scaling relations, (31)- 
(35) and (39)-(43), with which we can well understand what causes the differences among 
models. 

First, we have analyzed the dependence on the ambient number density. Our results 
are summarized as follows: 

1. At the typical ambient number density of GMCs, tih,o — 10 2 ~ 4 cm~ 3 , the molecular 
gas is finally accumulated in the shell, and gravitational fragmentation is expected. 

2. The timescale of the evolution, size of the H II region, and the swept-up mass in the 
shell quantitatively change with the different ambient number densities. With the lower 
number density, the H II region expands into the larger region over a long timescale, so 
that the swept-up mass in the shell is larger. The expected timescale within which the 
layer becomes gravitationally unstable is also longer with the lower ambient number 
density. 

3. The H2 molecules rapidly accumulate in the shell, but CO molecules are more sensitive 
to the shielding effect of FUV photons. With the lower number density, the accumu- 
lation of CO molecules in the shell is delayed owing to the small column density of the 
shell. At nn,o = 10 4 cm~ 3 , the shell is dominated by CO molecules before the frag- 
mentation. At nn,o = 10 2 cm -3 , however, the shell-fragmentation is expected to take 
place before the accumulation of CO molecules. If the shell fragments into clumps with 
high column density, CO molecules will be formed there owing to the self-shielding of 
clumps. 

Next, we have analyzed the dependence on the luminosity (mass) of the central stars. Our 
results are the following: 
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4. With the central stars of M* = 101 — 12 M Q , the molecular gas is finally accumulated 
in the shell, and gravitational fragmentation is expected. 

5. The timescale of the shell-fragmentation is relatively insensitive to the luminosity of 
the central star. At nn,o ~ 10 3 cm" 3 , the unstable region spreads in the shell over 
0.5—2 Myrs with the central star of M* = 101 — 12 M , whose UV and FUV luminosities 
differ by 2-4 orders of magnitude. 

6. The maximum density of the shell is always higher with the more massive central star. 
This is due to the difference in the Mach number of the SF. With the higher-mass star, 
the FUV radiation is easily shielded and the temperature in front of the SF quickly 
decreases, which raises the Mach number. 

According to the above results, we can conclude that an expanding H II region in homoge- 
neous molecular gas should be an efficient trigger for star formation. Once the massive star 
is ignited in the cloud, the star formation can propagate to the remaining molecular material 
around the exciting massive star. 

We are grateful to Shin Mineshige, Takashi Nakamura, Toru Tsuribe, Kazu Omukai, 
Hiroshi Koyama, and Akira Mizuta for their useful comments and encouragement. We also 
owe thanks to the YITP and NAOJ computer systems for the numerical calculations. SI 
is supported by the Grant-in-Aid (15740118, 16077202) from the Ministry of Education, 
Culture, Sports, Science, and Technology (MEXT) of Japan. 
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A. Finite Difference Method 

In this section, we explain the difference method for the ionization equation and the 
energy equation, (27). Other rate equations are differenced in the same way as the ionization 
equation. For simplicity, we temporarily omit all indices and drop the frequency- dependence 
of the radiation field in the ionization equation, 

*?L = (l-X)aF + R, (Ai) 

where R means some processes other than photoionization (e.g., recombination etc). The 
simplest difference form of (AI) is, 

XjW ~J = - {(1 - X 3 )F 3 T + ^ + (A2) 

where all physical quantities are defined at the cell center. The time average, Q/ +1/2 for a 
quantity, Q is defined as 

q; +1/2 = Q/+1 2 + Qr . (as) 

Using the radiative transfer equation, (8), the UV photon number flux in the j-th cell, Fj is 
written in time steps, 

r 2 , / r r i +i 



/•} /} L-^fexp y o j n(l-X)dr] (A I) 

Eliminating F- n+1 in equation (A2) with (A4), we can calculate X- n+1 iteratively with other 
rate equations and the energy equation. However, this differenced form provides no reliable 
results. If the optical depth through the j-th cell is thick enough, At, 1, the UV flux at 
the j-th cell center, Fj, significantly decreases from the UV flux at the previous cell center, 
Fj_i. Even if a sufficient number of UV photons reaches the cell interface at rj-1/2, only a 
few photons may reach the j-th cell center and the j-th cell will not be photoionized. This 
situation occurs at the ionization front (IF), where the optical depth changes sharply. 

The quick way to improve the scheme is to define the flux at the inner cell interface, not 
in the cell center. However, we adopt another more sophisticated differenced form for rate 
equations. With equation (Al) and radiative transfer equation (7), the ionization equation 
becomes, 

at nr A or 
First, equation (A5) is differenced for the spatial grids, 

dXi r 2 _ 1/2 F,_ 1/2 [l - exp(-n J -(r(l - X^Arj) 



dt w j(rJ+i/2 - r ?-i/ 2 )/ 3 



Rj ■ (A6) 
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In the right-hand side, the first term represents the consumption rate of UV photons for the 
photoionization per unit volume. This equation assures that the number of ionized atoms 
per unit time is equal to that of the consumed photons per unit time in the j-th cell, if other 
processes, Rj, are neglected. 

Since the ionization rate sharply changes at the IF, Xj should be calculated with some 
precision (e..g, Mathews 1965; Tenorio-Tagle 1976). We approximate the time evolution of 
the ionization rate between two time steps with the power-law form, 

j£ n+l X n 

X j (t)=X j n + j At j t , (0<t<At) (A7) 

and analytically integrate (A6) by time, t from t n to t n+1 . Other terms are calculated in the 
same manner. With the frequency-dependence of the radiation field and the dust extinction, 
we can derive the differenced form in a similar way. We use the same differenced form to 
calculate the photodissociation rate of H 2 and CO molecules. The self-shielding effect works 
with the low column density ; Nn 2 > 10 14 cm" 2 for H 2 molecules (Draine & Bertoldi 1996) 
and iVco > 10 15 cm -2 for CO molecules (Lee et al. 1996). The column density of one cell 
can be much higher than these values in our time-dependent calculations, and we do not use 
the shielding function given by Draine & Bertoldi (1996) and Lee et al. (1996). 

The energy equation (27) is differenced as, 



At \2w J \ At At 

+ 



Fh(7 ~ 1) f(u\ n+1/2 r 2 dp ( 1 V t+1/2 d(r 2 up) ( 1 V +1/2 d(r 2 q COQd ) 
k 1 \wJ dm \w ) dm \w ) dm 



7-1 /r-A 



k \ w 



n+l/2 



(A8) 



With differenced ionization equation (A6), other similar rate equations, and the energy 
equation (A8), we iteratively solve the new temperature, T- n+l and chemical compositions, 



A SJ - 



B. Necessity of Thermal Conduction — Field Condition — 

We include the heat transfer by thermal conduction in the energy equation (3). Recently, 
Koyama & Inutsuka (2004) have shown that the thermal conduction must be included to 
calculate the gas-dynamics with the thermal instability (TI). The conduction stabilizes the 
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TI for the short wavelength and the critical wavelength is called the Field length, 

'KT\ 1/2 
nA 



^=1— > (Bl) 



(Field 1965) where K is the conductivity and A is the cooling rate per hydrogen atom. 
Koyama & Inutsuka (2004) have shown that the Field length must be resolved with several 
spatial cells to converge the calculation, increasing the number of spatial grids. In our 
calculations, the TI occurs around T ~ several x 10 3 K mainly at the IF. If we do not include 
the thermal conduction, unphysical oscillation spreads from the IF and the calculation does 
not converge. Therefore, we need to include the thermal conduction in the energy equation. 

Since the Field length is generally much shorter than the characteristic length scale 
(e.g., Stromgren radius) in our calculations, we adopt another available method to resolve 
Field length. We first determine the width of the cell, Ar,- and calculate the conductivity, 
Kj so that the Field length is longer than the given width of the cell. The index j means 
the difference step in space. If we set Af ~ 3Arj, such conductivity is calculated as, 

Kj = ^ (Ar,) 2 . (B2) 

The conductivity, Kj is calculated every time step with equation (B2). If we set 100 cells per 
initial Stromgren radius for the ambient number density, 10 3 cm" 3 , the calculated conduc- 
tivity is, K ~ 4 x 10 11 " 10 erg cm" 1 K" 1 s -1 , which is significantly larger than the standard 
value. However, the conduction with equation (B2) smooths out the temperature structure 
only on the Field length scale, which covers a few cells (Af ~ 3 Arj). Since the meaningful 
temperature profile changes over a much longer length scale, we use the equation (B2) to 
calculate the conductivity. 



C. Doppler Shift of Lines 

C.l. Doppler Effect for Cooling Lines 

The dominant cooling processes in the PDR are the radiative loss by the fine-structure 
transitions, [OI] 63.1 /im and [CII] 157.7 //m (Hollenbach & Tielens 1999). These lines can 
be optically thick at iV H > 10 21 cm" 2 (Hollenbach & McKee 1979, 1989), and we use the 
escape probability for these lines (de Jong et al. 1980; Tielens & Hollenbach 1985). Since the 
swept-up shell expands at a supersonic velocity relative to the ambient gas, a few - 10 km/s, 
we take account of the Doppler shift for the trapping effect. The gas velocity in the shell 
is usually almost constant, and we separately calculate the escape probability for the gas in 
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the shell and ahead of the shell. For the gas in the shell, photons can escape from the inner 
(IF) and outer edge (SF) of the shell, 

Pesc — fesc( T x,sl) + /esc ( T x,s2) , (CI) 

where j3 esc is the escape probability, r XjS i and t x>s2 are the optical depth averaged over the 
line to the IF and the SF respectively, and / csc (t) is defined as, 

r<7, 

t>7 (C2) 

(de Jong et al. 1980). For the gas ahead of the shell, we assume that only half of the photons 
can escape from the SF, 

Pesc /esc v a:, ad) j (^3) 

where r x ^ is the optical depth to the SF. The CO molecules are another important coolant in 
the PDR. We do not solve the non-LTE rotational level population, to save the computational 
time, but adopt the analytical formulae given by Hollenbach & McKee (1979) including the 
trapping effect. We take account of the Doppler shift for the trapping effect in the same 
manner as the fine-structure lines. In our calculations, the rovibrational lines of H 2 molecules 
are reliably optically thin. 



/esc(r) 



1 - exp(-2.34r) 



4.68r 



4r 



In 



T 



~l 0.5 



C.2. Doppler Effect for Photodissociating Lines 

The FUV photons emitted from the central star photodissociate the molecules, both 
in the swept-up shell and in the unshocked ambient medium. However, we do not in- 
clude the Doppler effect for the photodissociating lines of H 2 and CO molecules. Be- 
low, we briefly verify this approximation. The Doppler shift in the FUV energy range 
is, hAu d = hu (u s \/c) ~ 10~ 4 eV, where hu ~ 10 eV is the characteristic energy of FUV 
photons, u s i ~ 10 km/s is the typical velocity of the shell, and c is the speed of light. The 
Doppler shift is larger than the Doppler width of lines owing to the supersonic velocity of 
the shell. As noted in §2.1, however, the absorption of FUV photons in the Doppler core is 
not important. In our calculations, for example, the column density of H2 molecules easily 
exceeds Nu 2 ~ 10 19 cm -2 , and the FUV photons are mainly absorbed in the Lorentz wings 
after that. Since the energy range of the Lorentz wings is much wider than that of the 
Doppler core, the Doppler shift of the core is not so important. 
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Fig. 1. — Our representative lines for Lyman-bands, 0Vh 2 - The points indicate the center of 
the grids. The dashed line means the dust absorption cross section for H 2 dissociating FUV 
photons, era 
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Fig. 2. — Our representative line for the photodissociation of CO molecules, Zqct^co- The 
points indicate the center of the grids. The dashed line (dot-dashed line) means the cross 
section for the dust absorption (shielding by H 2 molecules), o"d,co (°"h 2 ,co/2)- 
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Fig. 3. — The snapshots of the gas-dynamical evolution of model HD-S41. In each panel, five 
snapshots represent the profiles at t = 0.023, 0.069, 0.115, 0.161 and 0.207 Myr respectively 
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Fig. 4. — Same as Fig3 but for model LD-S41. In each panel, five snapshots represent the 
profiles at t = 0.83, 2.5, 4.2, 5.8 and 7.5 Myr respectively. 
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Fig. 5. — Upper panel: The time evolution of the various front positions in model HD-S41. 
The solid and dotted line mean the position of the IF and SF. The broken (dot-solid) line 
represents the DF of H 2 (CO), where X H2 = 2nn 2 /nH nuc = 0.5 and X C o = «co/nc mc — 0-5. 
We plot the position of DFs only before the SF catches up with each DFs. After the DFs 
are taken in the shell, the positions of DFs are close to that of the SF. Lower panel: The 
time evolution of the column density of each region in model HD-S41, where the vertical 
axis, Nr is the column density of the hydrogen nuclei included in all species of H + , H and H2 
from the central star. Thin contour lines represent the position where Xh 2 (dashed contours) 
and Xqo (dot-dashed contours) = 0.2, 0.5, and 0.8. The shaded region corresponds to the 
expected unstable region, where t > (Gp) _1//2 . 
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Fig. 6. — Same as Fig. 5 but for model LD-S41. 
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(1 ) High density case (2) Low density case 
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Fig. 7. — Schematic figures for the H II region and the swept-up shell in models with different 
ambient number densities. These figures show the snapshot at one t/£d y n in both models. 
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Fig. 8. — The density and the gas/dust temperature profiles in model HD-S41. These 
snapshots are at t — 0.161 Myr, which corresponds to ~ 15 tdyn- 
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Fig. 9. — Same as Fig. 8 but for model LD-S41. These snapshots are at t — 3.5 Myr, which 
corresponds to ~ 15 t& yIl . 



-42- 




Fig. 10. — Same as Fig3 but for model S101. In each panel, five snapshots represent the 
profiles at t — 0.1, 0.3, 0.5, 0.7 and 0.9 Myr respectively. 
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Fig. 11. — Same as Fig3 but for model S19. In each panel, five snapshots represent the 
profiles at t = 0.15, 0.45, 0.75, 1.05 and 1.35 Myr respectively. We add one early snapshot 
of the density profile at t = 0.02 Myr. 
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Fig. 12. — Same as Fig3 but for model S12. In each panel, five snapshots represent the 
profiles at t — 0.25, 0.75, 1.25, 1.75 and 2.25 Myr respectively. We add two early snapshots 
of the density profile at t = 0.025, and 0.125 Myr. 
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Fig. 13 



. — Same as Fig. 5 but for model S101. 
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Fig. 14. — Same as Fig.5 but for model S19. 
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Fig. 15. — Same as Fig.5 but for model S12. 
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Fig. 16. — Schematic figures for the H II region and the swept-up shell in models with 
different central stars. These figures show the snapshot at one t/t dyn in both models. 
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Fig. 17. — Upper panel : Time evolution of the maximum density of the shell in model S101 
(dotted line), S41 (solid line), S19 (dot-solid line) and S12 (broken line). Lower panel : Time 
evolution of the Mach number of the SF in each model. 



